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The quantum dynamics of a two-state system (qubit) can be governed by means of 
external control parameters present in time-dependent bias pulses of special forms. 
We consider the class of biases for which the time evolution equation without a 
dissipation can be solved exactly. Concentrating for deflniteness on the flux qubit we 
calculate the probability of the definite direction of the current in the loop and its 
time-averaged values as functions of the qubit 's control parameters both analytically 
and solving numerically the equation of motion for the density matrix in the presence 
of relaxation and decoherence. It is shown that there exist such time-dependent biases 
that the definite current direction probability with no dissipation taken into account 
becomes a monotonously growing function of time tending to a value which may 
exceed 1/2. We also calculate the probability to find the system in the excited state 
and show the possibility of the inverse population in a properly driven two-state 
system provided the relaxation and dephasing rates are small enough. 

PACS numbers: 03.67.Lx, 03.75.Lm, 85.25.Am 



I. INTRODUCTION 

In the past few years Josephson junctions-based devices have been widely studied 
both theoretically and experimentally as possible candidates for the implementation 
of a quantum computer I n fact, under appropriate values of external bias 

pulses, they behave as two-state systems which can be used as a model for quantum 
bits (qubits). Several systems like ion traps and NMR systems [n3], [ll] have been 
suggested for physical realizations of the qubit but Josephson devices being scalable 
up to large numbers of qubits as nanocomponents embedded into an electronic circuit 
exhibit the main technological advantage. Moreover, it is possible to prepare these 
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devices in a prefixed initial state or in a superposition of states and to control their 
dynamics by an external voltage and magnetic flux [12]. 

One common approach to control the qubit dynamics is to drive the two-state 
system, i.e. a particle in a double well potential, with an oscillating field. As a result 
an interesting phenomenon may take place. Instead of oscillating between the wells the 
particle may become localized in one of them. Such an unusual behavior of the quantum 
particle is known in the literature as coherent destruction of tunneling in double-well 
potential studies [l3]-[l8|, dynamic localization in transport analysis 19], 20] and 



population trapping in laser-atom physics 



2l|]. It is worth stressing that up to now 



the phenomenon was essentially related with an oscillating character of the driving 
external field. 

In the present article we show that the oscillating character of the field is not 
compulsory for appearance of the phenomenon. We present a class of non-periodical 
time-dependent bias pulses leading to a similar behavior of the qubit. 



A qubit in the two-state approximation can be described (see e.g. 
Schrodinger equation 



li, [221]) by the 



with the Hamiltonian 



H = Aa x + e (t) a z 



h 



(1) 



written down in the basis of "physical" states {|0) , |1)}, which are the eigenstates of the 
Pauli matrix a z (a z |0) = |0), er z |1) = — |1)) and \1/ = (ipi,ip2) T ■ In the case of a charge 
qubit |l|, these states correspond to a definite number of Cooper pairs on the island 
(Cooper-pair box). For a flux qubit [231 ] . they correspond to a definite direction of the 
current circulating in the ring. In this paper we assume that the tunnelling amplitude A 
is constant and the bias e is a time-dependent function, e = e(t). Bias e(t) is governed 
by gate voltage V g (t) of the gate electrode close to the Cooper-pair box in the case of a 
charge qubit and by magnetic flux $ x (t) piercing the qubit's loop in the case of a flux 
qubit. 

For definiteness we will consider here flux qubits. Then "physical" states {|0) , |1)} 
are the states with the definite (clockwise or counter-clockwise) direction of the current 
in the loop. The above mentioned time-dependent biases are, in fact, potentials for 
which Schrodinger equation (pp) can be solved exactly j^-j^zl]. Thus, the probability 
calculated using these exact solutions is, for example, probability of clockwise 
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current direction. We demonstrate that for some special non-periodical forms of 
potentials the probability of the clockwise current direction at the moment t, if at t = 
it was counter-clockwise, becomes a monotonously growing function of time tending 
to 3/4. Of course this is a strictly fixed excitation regime but we also study the 
behavior of the probability under small deviations from this specific regime. It should be 
noted that when parameters of the model are close enough to their specific values, the 
probability oscillates but its minimal value may exceed 1/2. It is established that such 
an unusual behavior of the probability is possible even in the presence of a dissipation. 
We study not only the time evolution of the probability, but also its time-averaged 
value. Moreover, we also calculate the probability of finding the system in the excited 
state and show the possibility of the inverse population in the two-state system even 
in the presence of a dissipation. Our main result is that using a properly chosen non- 
periodical time-dependent potentials one can "freeze" the qubit state i.e. localize only 
one of the two possible qubit states for a long time interval. We note that the probability 
of the definite current flow (or equivalently, the definite magnetic moment of the qubit) 
is related directly to experimentally measurable values such as the phase shift of the 



resonant circuit, weakly coupled to the qubit (as discussed, e.g., in Ref. (28||). Moreover, 
this "physical" basis is usually used in quantum computations. Thus, we hope that our 
results offer new opportunities for controlling the qubit behavior. Additional discussion 



about controlling the level population can be found in Refs. 



29|-(32j 



II. EXACTLY SOLVABLE BIAS PULSES 

In order to construct an exact solution of a differential equation the intertwining 
operator technique sometimes may be useful. The idea of the method dates from 



Darboux papers 



33| and is wide 



v used in the soliton theory Its quantum 



mechanical application (see e.g. [35|) is related to the fact the one-dimensional 
stationary Schrodinger equation is an ordinary second order differential equation 
defined by the operator of the potential energy. The method is based on the possibility 
to find an operator (intertwining operator) that relates solutions of two Schrodinger 
equations with different potentials. Thus, if one knows solutions of the Schrodinger 
equation with a given potential and an intertwining operator is available, there exists 
a possibility to construct solutions of the same equation with another potential, which 
cannot be completely arbitrary but is an internal characteristic of the method. There 



exists also a matrix-differential formulation of the method [34]] which was adapted to 



quantum mechanical problems in Ref. 361. 

n 

In [24j] it was shown how to construct differential-matrix intertwining operators 



for the system of two differential equations of type (pQ). For that authors [24] first 
reduce system (Dp) to the one-dimensional stationary Dirac equation with an effective 
non-Hermitian Hamiltonian where the time plays the role of a space variable and 
then apply the known procedure developed in [361 ] . Starting from the simplest case 
corresponding to s = eo = const, a new kind of nontrivial potentials (biases) for which 
Schrodinger equation (pp) can be solved exactly were found. Here we apply results 



obtained in [24|-[27j to describe the time evolution of the qubit, time dependence of 
the qubit localization probability and calculate its time-averaged value. 
Consider first the case when bias e = e(t) changes in the following way: 

In Ref. [24I a detailed analysis of solutions to equation (pQ) in this case is given. Therefore 
imposing the initial conditions (0) | 2 = 1 and \ip2 (0) | 2 = one finds probability 
(t) = \ip2 (t) | 2 of, for example, the clockwise current direction at the moment t if 
at t = it was counter-clockwise. For r = At and £ = ^ it reads 

P i (r) = 66(l + 4er 2 ) X 
"l6£ 4 6 2 T 2 cos 2 6r + 4£ 2 9r (l - 3£ 2 ) sin26r + (4£ 2 6 4 r 2 + (l - 3£ 2 ) 2 ) sin 2 6r , 

(3) 

It is clearly seen from here that P/ (r) is an oscillating function provided £ 2 7^ |. 
For £ 2 = I the probability becomes equal 

which is a function monotonically growing from zero at the initial time moment till the 
value 3/4 at r > 1 or ( > A -1 (see the thick line in Fig. [Th). It is important to note 
that at £ 2 close enough to | the value of probability P/ exceeds 1/2 very quickly (see 
the thin and dotted lines in Fig. (T^i). 

For the time-averaged probability one gets 
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It follows from Eq. (J5j) that at £ 2 = | the averaged probability exhibits a kind of the 
resonance behavior, i.e. it peaks to its maximal value P{ ~ 0.78, (see the thick line 
in Fig. [2]) although for any given £ P} (r) is a function asymptotically oscillating with 
frequency 4y| (see Fig. [IJ). 



More exactly solvable potentials may be obtained with the help of chains of the 
above simple transformations. Ref. j3| contains a detailed analysis of the properties of 
such chains. The authors show how to generate a large family of new exactly solvable 
biases for equation (DQ). For a two-fold transformation leading to the bias of the form 

e (45 - 180e 2 t 2 - 144e 4 t 4 + 64^t 6 ) 
£2[) 9 + 108egt 2 + 48^t 4 + 64^t 6 1 j 

the clockwise current direction probability is given by 

t 1 



10 (9 + 108£ 2 r 2 + 48£ 4 r 4 + 64f V 

x 



f 16£ 4 6 2 r 2 [I6£ 4 e 4 r 4 + 24£ 2 (3 - 14£ 2 + 7£ 4 ) r 2 + 9 (9 - 6£ 2 + £ 4 )] 
+ Qi [Q 2 sin 2 6r + Q 3 sin 20r] J (7) 

where 

Qi = 1 - m 2 + 5£ 4 , 

Q 2 = 64£ 6 6 4 r 6 + 48£ 4 (1 - 18£ 2 - 19£ 4 ) r 4 + 36£ 2 (3 - 2£ 2 + 11£ 4 ) r 2 + 9Qi , 
Q 3 = 12£ 2 6r (0 2 (16£ 4 r 4 + 7) + 2 (4£ 2 r 2 + 1) (1 - 5£ 2 )) . 

The last term in Eq. (j2|) describes the oscillations with frequency 26 and, hence, 
under the condition Qi — the oscillations in the time-dependence of the probability 
disappear and it acquires a monotonous character. Therefore, for e as given in (jHJ), 
unlike (j2J), we can indicate two possibilities for Q\ = 0. So, the probability of the 
clockwise current direction turns from an oscillating to a monotonous function of time 



both at £ = y 1 — 2/yfE and £ = y 1 + 2/ y/E. This behavior is illustrated in Fig. [TJd 
(thick lines) where we also show an oscillating character of the probability for the 
parameter £ close to the above critical values (thin and dotted lines). 
Time-averaged probability (|7j) for e (t) of the form ([6]) is given by 

— r 1 - 2£ 2 + 13£ 4 

p l = ir- • 8 

2 2 (1 + e 2 ) 

The ^-dependence of Pj is demonstrated in Fig.Eb by the thick line. It has a maximum 
P\ ss 0.91, at i « 1.46. 
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Let us consider a more complicated case j^J], [27] when the bias contains three free 
parameters 

e 3 (t) = eo + -r — , b 2 =e 2 -u 2 >0. (9) 

ocos (2ut + ip) — e 

Here <p is arbitrary but e and uo should satisfy the inequality given in Eq. ([9]). We 
note that in this case e = Ss(t) is a periodical function with an amplitude related with 
frequency. Expression ((9l) presents a generalization of formula ((2j). Indeed, putting 
ip = arctan ^ — | arctan j in Eq. j9]) , in the limit uj — > one recovers for £ result (j2J) • 
The analytic expression for -PJ(t) is rather involved and we will restrict ourselves 
to a graphical illustration of the clockwise current direction probability at ~ (see 
Fig. [3j). More graphical illustrations can be found in 



III. INFLUENCE OF A DISSIPATION ON PROBABILITIES 

A quantum system described by a wave function which is a solution of the 
Schrodinger equation with Hamiltonian ([TJ) interacts only with an external field 
described by function e(t). But in real experiments there is also an interaction of a 
Josephson junction device with an external reservoir which makes impossible describing 
the system in terms of a wave function since its state is not a pure quantum state 
anymore and should be described withe the help of a density matrix (operator in 
general see e. g [371 ]) formalism. Such kind of systems are particularly important in the 
context of quantum information processing where environment-induced decoherence is 
viewed as a fundamental obstacle for constructing a quantum information processors 
(e.g., fl). 

In this section we study the behavior of the qubit with bias e(t) as given in 
Eqs. ([2j), and taking into account the dephasing and the relaxation processes. We 
also study a possibility of the inverse population in the two-level system and investigate 
how it is influenced by a decoherence. It is worth noticing that in contrast to Ref. [291 ] . 
where the authors investigate a three-level system, we study a possibility of the inverse 
population for the two level system itself thus showing the possibility of building a 
two-level based laser working at a low temperature. 

To study the dynamic behaviour of the qubit, we use the master equation for the 
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density matrix. For the density matrix of the form 

we solve the equation of motion 

to obtain the probability P^ = 
on the system due to a weak coupling to the environment can be phenomenologically 
described by two parameters, the dephasing (T v ) and relaxation (T re i ax ) rates (see e.g. 
Ref. [39I]), which we introduce phenomenologically thus obtaining the following system 
of equations for X(t), Y{t) and Z(t): 

Z(0)). 

In order to verify the possibility of the inverse population in the two-level system we 
also calculate the probability P + [3^] of finding the system on the upper level (excited 
state). 

To see the influence of the relaxation on probabilities P/ and P+ exhibiting a 
monotonous time dependence for e(t) in form (121) we choose £ = l/Vo and plot them 
in Fig. [H We do not show the behavior of Pi and P 2 + corresponding to e(t) as given 
in since it is similar to that displayed in Fig. 0J Fig. El illustrates the evolution of 
probabilities P 3 T (t) and P 3 + when e{t) has form (jSJ) with £ = -^,/3 = ^p,^ = 0. Thick 
lines on these figures correspond to = T re i ax = when no relaxation is present in 
the system. Thin and dotted lines just show the relaxation effect for r^, = T re i ax = 0.01 
and = r re i ax = 0.1 respectively. All values are in units of A. It is clearly seen from 
here that the inverse population is still possible even if the dissipation is present but 
during a small time interval r only. 

Time-averaged probabilities P/ and Pj are plotted in Fig. Eh and [2}d as functions 
of dimensionless parameter £. Here the thick lines also correspond to the absence of 
the relaxation and the thin and dotted lines illustrate the relaxation effect for the 
same values of = T re i ax as in Figs. [4] and El The dash-dotted curves are drown for 
T^/A = T re i ax /A = 0.001. These results show the possibility of the inverse population 



l + Z X-iY 
X + iY 1 - Z 



dp 



[1 — Z(t)]/2. The effect of the relaxation processes 



^ = -2e(t)Y - r v X, 



-2AZ + 2e(t)X 



dY 
~dt 

— 21XY — T re i ax [Z 
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only for small enough dissipation (see the dash-dotted lines) and for parameter £ close 
enough to its critical values when the averaged probabilities pick to their maximal 
values. 

IV. CONCLUSION 

In conclusion, a two-state system subjected to biases of special forms, for which the 
Schrodinger equation is exactly solvable, was considered. For definiteness we studied 
the time-evolution of the Josephson-j unction qubit, but the similar consideration can 
be applied to other realizations of a two-state system. We have demonstrated that 
the dynamics of the Josephson-j unction qubit displays in these cases several nontrivial 
features. Varying the form of the time dependent bias, one can change qualitatively 
the dynamic behavior of the occupation probabilities. In particular, the amplitude of 
the oscillating probability, describing the definite current direction in the loop, can 
be tuned to zero, thus turning the probability to a monotonous function of time. We 
calculate the probability of finding the system in the excited state and show that 
the inverse population in the two-state system is possible even in the presence of a 
dissipation. The observation of such an occupation probability behavior may be related 
to experimentally measurable values, which would make an experimental verification 
of the theoretical predictions of the present paper. 
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Phc. 1: Time dependence of clockwise current direction probabilities, (a) P{ probability 
evolution at £ = iA (thin line), £ = 1 (dotted line) and £ = w | (thick line), (b) 



probability evolution at £ = w2 + £ = wl.2 + -J= and £ = wl + on the upper plot, 



£ = y/l.05 — -^g, £ = wl.01 — ^= and £ = wl — on the lower plot (dotted, thin and thick 
lines respectively). 




Phc. 3: Time dependence of clockwise current direction probabilities, ig probability evolution 
at 6 = 6.88, 9 = 7, £ = v 7 ^, = 0. 
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Phc. 4: The influence of relaxation on the probabilities, (a) The clockwise current direction 
probability P} at £ = (b) The upper level occupation probability P±~ at £ = l/v3^ 
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Phc. 5: The influence of relaxation on the probabilities, (a) The clockwise current direction 
probability P^ at £ = -^=, b = -^p, tp = 0. (b) The upper level occupation probability Pg~ at 
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